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Methods 

A set of speckle-tracking algorithms"' was used to determine the 1992, 1994, 1995 and 2000 
velocities from 1-24-day image pairs. Speckle tracking uses the displacements of the 
correlated speckle patterns in pairs of SAR images to derive ice motion estimates. 
Individual errors were up to a few hundred metres per year {see Fig. 2), but errors on 
averages {for example, Fig. 3) are below 100 myr . We did not tide-correct the speckle- 
tracked data, so there are biases on the floating ice that do not spatially average out. To 
assess this error, we estimated velocity for five 1992 InSAR pairs, each with different tidal 
errors. The standard deviation for these estimates was 69 myr . Our 1992 and 1994 
estimates are temporal averages of multiple (2 to 5) same-year pairs, which further reduces 
this error. The 2001 through 2003 estimates were derived using the IMCORR 25 feature- 
tracking software applied to 16-to-64-day Landsat image pairs. Established methods 26 
were applied to passive microwave data to determine the 2002 melt extent. 
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The summer of 2003 was probably the hottest in Europe since at 
latest ad 1500 1 " 4 , and unusually large numbers of heat-related 
deaths were reported in France, Germany and Italy 5 . It is an ill- 
posed question whether the 2003 heatwave was caused, in a simple 
deterministic sense, by a modification of the external influences on 
climate — for example, increasing concentrations of greenhouse 
gases in the atmosphere — because almost any such weather event 
might have occurred by chance in an unmodified climate. However, 
it is possible to estimate by how much human activities may have 
increased the risk of the occurrence of such a heatwave 6 " 8 . Here we 
use this conceptual framework to estimate the contribution of 
human-induced increases in atmospheric concentrations of green- 
house gases and other pollutants to the risk of the occurrence of 
unusually high mean summer temperatures throughout a large 
region of continental Europe. Using a threshold for mean summer 
temperature that was exceeded in 2003, but in no other year since 
the start of the instrumental record in 1851, we estimate it is very 
likely (confidence level >90%) 9 that human influence has at least 
doubled the risk of a heatwave exceeding this threshold magnitude. 

Temperatures near the Earth's surface are rising globally 10 , and 
evidence is mounting that most of the warming observed in recent 
decades has been caused by increasing atmospheric concentrations 
of greenhouse gases 9 ' 11 ' 12 . Anthropogenic increases in annual-mean 
temperatures have also been detected on continental scales, in 
Europe, North America and other land regions 13 " 15 . We first inves- 
tigate the origins of long-term changes in decadal-mean European 
summer (June-August) temperatures, determining the changes 
attributable to anthropogenic drivers of the climate system and 
changes attributable to natural drivers. We then estimate how the 
risk of mean June-August temperatures exceeding a particular 
extreme threshold in any individual summer has changed as a result 
of this anthropogenic interference in the climate system. 

Over the course of the twentieth century, June-August tempera- 
tures in Europe exhibited an overall increase, and a distinctive 
temporal pattern of temperature change, including cooling in the 
1950s and 1960s (Fig. 1). We focus on the region bounded by 10° W 
and 40° E and 30-50° N (Fig. 1 inset), this being one of the regions 
chosen in previous studies 1316 to represent climatically coherent 
regions sufficiently large to exhibit climate change signals above the 
noise of natural internal variability. We use a pre-selected region in 
order to minimize any bias that could result from selecting our 
region already knowing where the most extreme temperatures 
occurred. Even in such a large domain, 2003 was the warmest 
summer on record. The history of temperature change averaged 
over this region is well reproduced by simulations of the HadCM3 
climate model 17 , even at the model's relatively low spatial resolution 
(3.75° longitude by 2.5° latitude), when driven with both anthropo- 
genic and natural drivers of climate change (Fig. 1; see red, green, 
blue and turquoise lines). Four simulations (denoted ALL) were 
made with different initial conditions 18 , each with the same com- 
bination of well mixed greenhouse gases, sulphate aerosols and 
changes in tropospheric and stratospheric ozone, as well as natural 
changes in solar output and explosive volcanic eruptions 12 . A 
calculation of the temperature changes due to natural drivers 
alone (obtained by combining a simulation with solar forcing and 
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Figure 1 June-August temperature anomalies (relative to 1 961-90 mean, in K) over the 
region shown in inset. Shown are observed temperatures (black line, with low-pass- 
filtered temperatures as heavy black line), modelled temperatures from four HadCM3 
simulations including both anthropogenic and natural forcings to 2000 (red, green, blue 
and turquoise lines), and estimated HadCM3 response to purely natural natural forcings 



(yellow line). The observed 2003 temperature is shown as a star. Also shown (red, green 
and blue lines) are three simulations (initialized in 1 989) including changes in greenhouse 
gas and sulphur emissions according to the SRES A2 scenario to 21 OO 22 . The inset shows 
observed summer 2003 temperature anomalies, in K. 



a simulation with volcanic forcing, and denoted NAT) shows no 
warming in the latter part of the century (Fig. 1, yellow line). 

A necessary requirement for any detection of significant warming 
is that HadCM3 adequately represents the natural internal varia- 
bility of European summer temperatures. A quantitative test, in 
which an estimate of the forced response calculated from the ALL 
ensemble mean is subtracted from the observations and the residual 
compared with HadQVB simulated internal variability, shows no 
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Figure 2 Power spectra of European mean summer temperature. Solid line, observed 
spectrum after removing an independent estimate of the externally forced response 
provided by the ensemble mean of the ALL simulations. Shaded region, 5 to 95 percentile 
region of the estimated range of spectra of natural internal variability estimated from 
segments taken from the HadCM3 control run of the same length as the observations. 
Model spectral densities have been inflated by a factor of 1 .25 to allow for sampling 
uncertainty in the ALL ensemble mean. 
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significant discrepancy in variance on either interdecadal or 2- to 
10-year timescales (Fig. 2). If anything, the model appears to be 
overestimating observed variability in this season and region, 
making our estimates of attributable risk relatively conservative. 
We also find no evidence of a secular change in variance on sub- 
decadal timescales in either model or observations: over the twenty- 
first century, the standard deviation of the forced simulations under 
the SRES A2 scenario (Fig. 1), when a second-order polynomial 
trend has been removed, increases by an insignificant 0.01 K from 
the first 40 years (1990-2030) to the last 40 years (2060-2100). This 
contrasts with the conclusions of ref. 2, although their results 
applied to a much smaller region. 

We now apply a standard optimal detection analysis to European 
summer temperatures, similar to those applied to global scale 
patterns of temperature change 1112,19 . The analysis is a regression 
between decadal-mean seasonal-mean observations over 1920-99 
and simulated temperature changes over the same period from 
both anthropogenic and natural forcings (ALL) and from natural 
forcings alone (NAT), and is in most respects identical to that of 
ref. 13 (see Methods). Figure 3a shows the estimated likelihood 
functions for the factors by which we could scale up or down the 
amplitudes of the model-simulated responses to anthropogenic 
forcing (red curve) and natural forcing (green curve) while remain- 
ing consistent with the observations. As the fifth percentile of the 
scaling factor on the anthropogenic response (red curve in Fig. 3a) is 
greater than zero, an anthropogenic influence on decadal-mean 
European summer temperature is detected at the 5% significance 
level (that is, the hypothesis that there is no positive anthropogenic 
influence can be rejected at the 5% level). We conclude from this 
investigation of decadal-mean summer temperatures that it is very 
likely that past anthropogenic forcing is responsible for a significant 
fraction of the observed European summer warming. 

We now calculate the changed risk of extremely hot summers that 
can be attributed to past anthropogenic emissions, allowing for 
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uncertainty in the anthropogenic warming. Averaged over the 
region of interest (Fig. 1 inset), summer temperatures in 2003 
exceeded the 1961-90 mean by 2.3 K (Fig. 1, black star). To quantify 
changes in risk, we need an objective definition of the event in 
question. Using 2.3 K itself is problematic for three reasons: first, 
relying too closely on the details of what actually occurred when 
defining what we are looking for introduces a selection bias in our 
attribution procedure; second, temperature anomalies in 2003 may 
have been amplified by soil-moisture feedbacks 2 or interactions 
with the North Atlantic 20 , both of which may be under-estimated, 
although we do observe similar magnitude spikes in model summer 
temperatures in Fig. 1; third, given the length of model-simulated 
variability we have available, inferring the statistics of temperature 
excursions over 2 K requires extrapolation of extreme value distri- 
butions, which introduces further uncertainties. 2003 was the first 
year to reach or exceed a threshold of 1.6 K (2001 being the second- 
warmest European summer, at 1.5 K). We therefore consider 
how the probability of exceeding this threshold has changed, by 
comparing this estimated late-twentieth-century probability with 
the estimated probability of exceeding the same absolute threshold 
if there had been no anthropogenic influence on climate. Increasing 
this threshold to any value up to 2.3 K strengthens our conclusions 
regarding attributable risk; hence using a threshold that only just 
exceeds the second warmest summer is relatively conservative. 

Assuming that sub-decadal continental-scale variability is station- 
ary and adequately represented by HadCM3, we estimate possible 
distributions of temperatures in individual summers in the presence 
and absence of anthropogenic influence. We do this by adding 
HadCM3 control variability to reconstructions of 1990s decadal- 
mean temperatures both with all external factors included, and with 
anthropogenic factors removed, allowing for uncertainty in these 



decadal-mean temperatures from the detection analysis (Fig. 3b). 
Figure 4a shows the estimated likelihood of the risk (probability) of 
exceeding a 1.6 K threshold in the presence of anthropogenic climate 
change (red line) and in the absence of anthropogenic change (green 
line), expressed both as a frequency (number of occurrences per 
thousand years, bottom axis) and as a return period (top axis). The 
clear shift from the green to the red distribution implies that an 
appreciable fraction of the risk of such hot summers can be attributed 
to human influence on climate. Even in the presence of anthropogenic 
warming, we conclude that the estimated probability of exceeding 
1.6 K appears to be low (best estimate is a 1 in 250 year event (Fig. 4a, 
red curve) but this risk may be increasing rapidly). 

The fraction attributable risk (FAR) is estimated in Fig. 4b (see 
Methods). In certain circumstances, the figure of relevance in estab- 
lishing possible liability for compensation has been FAR = 0.5, 
corresponding to a doubling of risk over natural conditions 21 (mean- 
ing that one event in two would have happened naturally). According 
to our calculation, there is a greater than 90% chance that over half the 
risk of European summer temperatures exceeding a threshold of 1 .6 K 
is attributable to human influence on climate. Although there is a 
large spread, reflecting the remaining uncertainties in the effects of 
climate change on this spatial scale, the anthropogenic FAR could be 
substantially greater than 0.5. Also marked on Fig. 4b is a vertical line 
representing an overall 'best estimate' of the human contribution to 
the increased risk of these very hot European summers 7 , given the 
information that we have available at present. On this basis, human 
influence is to blame for 75% of the increased risk of such a heatwave. 

Our analysis shows that European summers are warming owing 
to anthropogenic climate change. Under un-mitigated emissions 
scenarios, summers like 2003 are likely to be experienced more 
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Figure 3 Estimated likelihood functions for anthropogenic and natural contributions to 
European summer temperature changes. The curves show estimated distributions of 
anthropogenic (red) and natural (green) scaling factors on model-simulated responses (a). 
1990s summer temperatures (relative to pre-industrial climate) including all external 
drivers of climate change (red) and with anthropogenic drivers removed (green) (b). A 
scaling factor of zero (horizontal solid line in a) implies no contribution to observed 1 990s 
temperatures from this driver, while unity (horizontal dashed line in a) implies no 
systematic under- or over-estimate by the model of the observed response to this driver. 
The width of these distributions reflects the uncertainties for these probabilities. 
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Figure 4 Change in risk of mean European summer temperatures exceeding the 1 .6 K 
threshold, a, Histograms of instantaneous return periods under late-twentieth-century 
conditions in the absence of anthropogenic climate change (green line) and with 
anthropogenic climate change (red line), b, Fraction attributable risk (FAR). Also shown, as 
the vertical line, is the 'best estimate' FAR, the mean risk attributable to anthropogenic 
factors averaged over the distribution. 
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frequently in future; HadCM3 projections (Fig. 1) indicate that the 
probability of European mean summer temperatures exceeding 
those of 2003 increases rapidly under the SRES A2 scenario 22 , 
with more than half of years warmer than 2003 by the 2040s. By 
the end of this century, Fig. 1 shows that 2003 would be classed as an 
anomalously cold summer relative to the new climate, for the 
scenario and model under consideration. 

We may have underestimated the FAR of a heatwave in 2003 
by using 1990s figures for attributable background warming. Con- 
versely, the FAR may have been overestimated by selecting an area 
that has recently been subject to extreme temperatures, although 
this effect should be alleviated by our use of an independently 
specified region. Modelling and forcing uncertainties, and any 
errors in the variability characteristics of the model, could mean 
that our assessment of the likelihood of the FAR exceeding 0.5 is in 
error. Including a different combination of the available natural 
HadCM3 simulations changes the results quoted by less than 
5%, but a systematic exploration of both modelling and forcing 
uncertainty would require a very large multi-model ensemble, 
varying model parameters across their range of possible values 
and exploring the range of potential forcing estimates 23 ' 24 . 

Anthropogenic warming trends in Europe imply an increased 
probability of very hot summers. Given the effects of the 2003 
heatwave, this suggests a greater risk of associated adverse impacts, 
although to make a quantitative attribution assessment of any 
specific social, economic or ecological impact will require detailed 
modelling of both local meteorological conditions and their 
relationships with the impact in question. Nevertheless, it seems 
likely that past human influence has more than doubled the risk of 
European mean summer temperatures as hot as 2003, and with the 
likelihood of such events projected to increase 100-fold over the 
next four decades, it is difficult to avoid the conclusion that 
potentially dangerous anthropogenic interference in the climate 
system is already underway. The FAR provides a potential measure 
for redistributing the associated costs of such extreme events 21 . □ 

Methods 

Attribution of decadal-mean seasonal- mean changes 

In a natural extension of previous work 12,13,23 , observed decadal-mean summer near- 
surface temperature changes, y, over land for 1920-99 are constructed from monthly- 
mean values extracted from the CRUTEM2(v) data set 10 for the region 10° W-40° E, 
30-50° N, requiring data for at least two out of three months { June-August) in at least half 
the years. They are regressed against the mean response over the observational data-mask 
of a four-member ensemble of HadCM3 simulations driven with both anthropogenic and 
natural forcings (xj, ALL, red, green, blue and turquoise lines in Fig. 1 } and a second set of 
simulations of the response to solar and volcanic forcing alone (x 2 , NAT, yellow line in 
Fig. 1) plus noise: 

y = (Xl - + (x 2 " V2 >& + v o (1) 

where (3 2 are the unknown scaling factors to be estimated in the regression, v Q is the 
noise in the observations, and c : and v 2 account for sampling uncertainty introduced by 
estimating model-simulated responses from a finite ensemble 26 . A scaling factor of 1 
would imply that the response in the model simulations is identical to the observed 
changes, while a factor of implies no correspondence between modelled and observed 
responses. Likelihood functions for scaling factors obtained from the regression account 
only for sampling variability that arises from natural internal climate variability as 
represented by HadCM3. They do not account for systematic errors in the shapes of the 
modelled patterns or uncertainties due to missing forcings (for example, land use changes 
or fossil-fuel black carbon emissions}. If we assume a uniform prior in these scaling 
factors, these likelihoods can be thought of as probability distributions. 

The combined effects of solar and volcanic forcings are estimated from a simulation 
with volcanic forcing amplified by a factor of 10 and solar forcing amplified by a factor of 5, 
multiplied by 0.1 and 0.2, respectively, and added. This approach was taken in order to 
obtain the clearest possible climate signal from a weak forcing and a limited number of 
climate simulations. Previous work showed that the large-scale temperature responds 
linearly to increasing these forcings 27 . Including an additional four simulations with 
unamplificd natural forcings in the analysis makes little difference to the results. 

We use the HadCM3 control run and intra-ensemblc variability to estimate internal 
variability v , v \ and v 2 , taking into account the enhanced signal-to-noise from averaging 
over the ALL ensemble and inflating the solar and volcanic forcing 13 . This gives scaling 
factors for the combined effects of all forcings and for natural forcings, from which the 
anthropogenic and natural scaling factors can be calculated by linear transformation 27 
(Fig. 3a). Where the 5 to 95 percentile uncertainty range does not include zero scaling 
factor, this indicates that the relevant signal has been detected at the 5% significance level, 
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that is, there is less than 5% chance that the relevant forcing has had no effect. From the 
scaling factors and their uncertainties, we derive likelihood functions of 1990s 
temperature anomalies attributable to the combined effects of all forcings and attributable 
to natural factors alone (Fig. 3b). 



Attribution of change in risk of exceeding a threshold 

The distribution of attributed anthropogenic decadal-mean warming in the 1990s is used 
to calculate the FAR for an individual European summer exceeding the 1.6 K threshold as 
follows. We estimate the probability of exceeding the 1.6 K threshold from the generalized 
Pareto distribution (GPD) fitted to the HadCM3 control run temperatures 28 . The spectral 
analysis shown in Fig. 2 indicates that HadCM3 provides an adequate representation of 
internal variability. We calculate the probability, P , of the 1.6 K threshold being exceeded 
without anthropogenic climate change (shown by the green curve in Fig. 4a, expressed as 
return period, top axis, and instantaneous risk: rate of occurrence per 10 yr, bottom axis) 
and the probability, P u of exceeding the threshold with anthropogenic climate change (red 
curve in Fig. 4a). In estimating the green curve, the chance of exceeding the 1.6 K threshold 
is re-estimated using the GPD with the mean summer temperature adjusted to each 
percentile of the estimated distribution of decadal-mean temperature anomalies 
attributable to natural forcing and internal variability. For the red curve, the response to 
anthropogenic forcing is also included. A bootstrap resampling method is used to allow 
for uncertainty in GPD parameters. From these distributions, we calculate the distribution 
of attributable increase in risk (iVP ) and the FAR ((P l — Po)/?^, as shown in Fig. 4b, 
whose spread reflects uncertainty in the magnitude of the forced decadal-mean response 
and in the estimation of the extreme value distribution parameters due to limited sample 
size. 
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Multidisciplinary investigations at the Los Ajos archaeological 
mound complex in the wetlands of southeastern Uruguay chal- 
lenge the traditional view that the La Plata basin was inhabited by 
simple groups of hunters and gatherers for much of the pre- 
Hispanic era 1 " 4 . Here we report new archaeological, palaeo- 
ecological and botanical data indicating that during an increas- 
ingly drier mid-Holocene, at around 4,190 radiocarbon ( 14 C) 
years before present (bp), Los Ajos became a permanent circular 
plaza village, and its inhabitants adopted the earliest cultivars 
known in southern South America. The architectural plan of 
Los Ajos during the following Ceramic Mound Period (around 
3,000-500 14 C yrBP) is similar to, but earlier than, settlement 
patterns demonstrated in Amazonia 5 " 10 , revealing a new and 
independent architectural tradition for South America. 

Research on the emergence of complex societies in South 
America has mainly concentrated on Andean coastal and highland 
valleys 11 " 14 , and more recently in the lowland forest and riverine 
regions of Amazonia 3 " 10 . The La Plata basin (Fig. la) is a large and 
little explored river system that is beginning to reveal an early and 
long sequence of unique and complex cultural trajectories. The 
natural environment of the study area is dominated by subtropical 
grasslands interspersed with vast extensions of wetlands. In strategic 
locations circumscribed to wetland floodplains, archaeological 
mound complexes are large and numerous. They have mounded 
architecture geometrically arranged in circular, elliptical and horse- 
shoe formats that surround a central communal space (Fig. lb). 

The first excavations at the multi-mound site called Los Ajos, by 
Bracco 15 , in the early 1990s consisted of a block excavation in 
Mound Alfa, a test unit in Mound Beta and a few opportunistic 
test units in off-mound areas. This work established the mid- 
Holocene age of the earthen mounds in the area. The Preceramic 
Mound Period (PMP) component at Los Ajos yielded five dates 
between 3,950 and 3,350 14 C yr bp (4,580 and 3,380 calibrated (cal.) 
yrBP) 16 . Closely comparable dates from the deeper PMP com- 
ponents of the Puntas de San Luis, Isla Larga and Potrerillo sites 
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collectively ascertained the antiquity of the PMP 17 " 19 (Table 1). 

Our new excavation programme consisted of the placement of a 
block excavation in Mound Gamma, a test unit in Mound Delta, 
two trench transects articulating mound and off-mound areas and a 
50-m systematic interval transect sampling strategy of test units to 
target off-mound areas (Fig. 2) totalling an excavated area of 305 m 2 
(ref. 20). Our work revealed that Los Ajos is one of the largest and 
most formally laid out sites in the study area and covers about 12 ha 
(Fig. 2a). Its Inner Precinct includes six flat-topped, quadrangular 
platform mounds (called 6, Alfa, Delta, Gamma, 4 and 7) closely 
arranged in a horseshoe formation and with a height above ground 
level of 1.75 to 2.5 m (Fig. 2a, b). Two dome-shaped mounds (called 
Beta and 8) frame the central, oval plaza with a size of 75 X 50 m 
(Supplementary Fig. la). The formal and compact Inner Precinct 
contrasts with more dispersed and informally arranged peripheral 
sectors, which include two crescent-shaped rises (named TBN 
(Supplementary Fig. lb) and TBS), five circular and three elongated 
lower dome-shaped mounds, borrow pits and a vast off-mound area 
bearing subsurface occupational refuse. 

A series of major social and economic changes took place at Los 
Ajos during the PMP. The broad contemporaneity of radiocarbon 
dates (Table 1), artefact content and similarities in Preceramic 
Mound Component (PMC) stratigraphy among mounds Alfa, 
Delta and Gamma suggest that the Los Ajos inhabitants began to 
live in a circular household-based community, partitioning the site 
into a number of discrete functional areas characterized by the 
placement of residential units around a central plaza area. Charcoal 
from the basal level of Mound Gamma, 270-275 cm deep (arbitrary 
depth), dates the beginning of the PMC at 4,190 ± 40 14 C yrBP 
(4,840-4,580 cal. yrBP). Another radiocarbon assay at 205-210 cm 
deep yielded a date of around 3,460 14 C yrBP (3,980-3,470 cal. 
yr bp). The upper portion of the PMC in Mound Delta yielded a date 
of around 2,960 ± 120 14 C yrBP. Taken together, the eight dates 
from Los Ajos place the PMC occupation between 4,190 ± 40 and 
2,960 ± 40 14 C yr bp. The two oldest dates from the basal levels of 
the PMC at Mound Gamma and Alfa suggest that mound-building 
began between around 4,190 and 3,950 C yrBP (4,840-4,160 cal. 
yr bp). Mounds grew as a result of multiple overlapping of domestic 
occupations where a wide range of activities associated with food 
preparation, consumption, stone tool production and maintenance 
took place. The PMC Layer 4 is characterized by a ~85-cm-thick 
compact, very dark brown silty loam organic sediment (Fig. 3) 
consisting of relatively undifferentiated deposits composed of lithic 
debitage and tools, small fragments of charred bone, ash and soot 
lenses, and small pieces of burned clay. 

Our associated palaeoecological data 20 indicate that, similar to 
other regions in the Americas 21 , the mid-Holocene (between 
6,620 ± 40 14 C yrBP (7,580-7,440 cal. yrBP) and 4,020 ± 40 14 C 
yrBP (4,570-4,410 cal. yrBP)) was a period of significant environ- 
mental flux marked by increasing aridity. At around 4,020 14 C yr bp 
(4,570-4,410 cal. yrBP) a maximum drying episode occurred, as 
evidenced by a massive spike of Amaranthaceae/Chenopodiaceae 
coupled with a sharp drop in wetland species (Supplementary 
Fig. 2). These changes probably caused a decrease in the surface 
water recharge to the inland wetlands and waterways. Although 
reduced in extent, wetlands probably became attractive places for 
pre-Hispanic populations by providing abundant, now more highly 
circumscribed plant and animal resources and a stable source of 
water. The mid-Holocene drying trend may thus have acted as an 
important catalyst for the reorientation of settlements towards 
the topographically higher freshwater wetlands, where permanent 
communities were established. 

Despite the application of an intensive flotation program, seeds, 
roots and nuts were not recovered; however, phytoliths and starch 
grains were abundant (Methods are available as Supplementary 
Information). Phytoliths diagnostic of maize cobs 22 (Fig. 4c) record 
their earliest appearance in the level from 255-260 cm deep, 15 cm 
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